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Abstract 

Flight 10 of NASA's Subsonic Assessment (SASS) Ozone and Nitrogen Oxide 
Experiment (SONEX) extended southwest of Lajes, Azores. A variety of chemical 
signatures were encountered. These signatures are examined in detail, relating them to 
meteorological data from a high resolution numerical model having horizontal grid spacing 
of 30 and 90 km and 26 vertical levels. The meteorological output at hourly intervals is used 
to create backward trajectories from the locations of the chemical signatures. 

Four major categories of chemical signatures are discussed — stratospheric, 
lightning, continental pollution, and a transition layer. The strong stratospheric signal is 
encountered just south of the Azores in a region of depressed tropopause height. Three 
chemical signatures at different altitudes in the upper troposphere are attributed to ligh tnin g 
Backward trajectories arriving at locations of these signatures are related to locations of 
cloud-to-ground lightning. Results show that the trajectories pass through regions of 
lightning 1-2 days earlier over the eastern Gulf of Mexico and off the southeast coast of the 
United States. The lowest leg of the flight exhibits a chemical signature consistent with 
continental pollution. Trajectories arriving at this signature are found to pass over the highly 
populated Northeast Corridor of the United States. Surface based pollution apparently is 
lofted to the altitudes of the trajectories by convective clouds along the East Coast that did 
not contain lightning. Finally, a chemical transition layer is described. Its chemical 
signature is intermediate to those of lightning and continental pollution. Trajectories arriving 
in this layer pass between the trajectories of the lightning and pollution signatures. Thus, 
they probably are impacted by both sources. 
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1. Introduction 

NASA's Subsonic Assessment (SASS) Ozone and Nitrogen Oxide Experiment 
(SONEX) was conducted over the North Atlantic Flight Corridor (NAFC) during Fall 1997 
to study the effects of subsonic aircraft emissions on the atmospheric chemistry of the upper 
troposphere and lower stratosphere. During the sixteen flights comprising the SONEX field 
phase, NASA's instrumented DC-8 aircraft collected in situ high resolution chemical data 
throughout the NAFC and surrounding areas. Singh et al. [this issue] give a comprehensive 
overview of the SONEX program. 

SONEX Flight 10 on October 29, 1997 was a tropical survey originating at Lajes, 
Terceira Island, Azores (38.7°N, 27.1°W). The flight track was oriented northeast to 
southwest, reaching a southwestern terminus near 19.9°N, 37.1°W (Figure 1). A stack 
maneuver of four legs was performed between the southernmost latitude and ~25°N to 
sample chemical characteristics at multiple levels. Although the goal of Flight 10 was to 
sample clean, subtropical air exhibiting minimal continental, convective, and aircraft 
influence, complex chemical signatures from a variety of sources actually were encountered. 

Insight into the transport history and origins of air parcels can be gained by relating 
meteorological data to chemical data. Species-species correlations are useful since they can 
indicate the types of atmospheric processes that air samples have experienced, e.g., 
convection, pollution, mixing, and photochemical reactions [Davis et al., this issue, 
Thompson et al., this issue a]. Trajectory calculations can indicate possible source regions, 
transport paths, and age of the air [e.g., Board et al., 1999]. 

Origins and transport paths of atmospheric chemicals over the North Atlantic Ocean 
are not easily investigated. This is due to the sparsity of meteorological data over oceanic 
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regions, as well as the vast distances that air parcels can travel over short periods of time in 
the mid-latitudes. Transport paths and consequently, the chemical composition of air in this 
region can vary greatly, particularly during the cool season due to strong winds and high 
amplitude flow patterns. 

This research relates observed chemical signatures during SONEX Flight 10 to high- 
resolution numerically derived meteorological data. Our objective is to understand the 
chemical signatures by determining source regions and transport paths of the air parcels that 
were sampled. We use the National Center for Atmospheric Research (NCAR)/Penn State 
Fifth Generation Mesoscale Model (MM5) [Grell et al . , 1994] to create a high resolution 
three dimensional set of meteorological data. Backward trajectories arriving at locations of 
the chemical signatures then are calculated using the model-derived winds. Finally, the 
trajectory paths are related to meteorological and/or anthropogenic chemical sources. 

We believe that our high resolution MM5-derived meteorological data (30 or 90 km, 
hourly) provide a better representation of actual wind regimes and, consequently, more 
accurate trajectories than is possible using global gridded analysis (typically 1 10 km, 6- 
hourly). This improved resolution is important since trajectories are highly sensitive to small 
errors in both wind direction and speed as well as the temporal frequency of the wind data 
[Doty and Perkey, 1992]. Other studies that have used high resolution numerically-derived 
meteorological data to compliment the synthesis of high resolution chemical data include 
Cha field et al. [1996, this issue], Pickering et al. [1996, this issue], Liu et al. [1996], and 
Wangetal. [1996]. 
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2. Data and Methodology 
2.1 Meteorological Data 

We employed a two-way interactive nested configuration for the non-hydrostatic 
MM5 simulations (Figure 2). The coarse outer grid extended far beyond the SONEX region, 
covering most of North America and Europe, to reduce the effects of boundary errors 
[Warner et al., 1997]. The horizontal spacing of the coarse grid was 90 km, with 147 grid 
points in the east to west direction and 87 grid points south to north. The higher resolution 
inner grid, centered over the path of Flight 10, had a horizontal spacing of 30 km. This gnd 
contained 136 points east to west and 166 grid points south to north. Both domains utilized 
26 sigma levels in the vertical that were relatively closely spaced near the tropopause (-20 
hPa separation). The Grell cumulus scheme [Grell, 1993] and the Blackadar planetary 
boundary layer scheme [Zhang and Anthes, 1 982] were chosen as the model physics 
parameterizations. 

The MM5 was initialized with a gridded data set prepared by the European Centre for 
Medium-Range Weather Forecasting (ECMWF) [Bengtsson, 1985, Hollingsworth et al., 
1986]. The resolution of these data was 2.5° by 2.5° in the horizontal with 15 isobaric levels 
in the vertical between 1000 and 10 hPa. Boundary conditions at subsequent times also were 
obtained from the ECMWF data. Sea surface temperatures were taken from a data set 
provided by the National Center for Environmental Prediction (NCEP). The 67 hour 
simulation on the 90 km domain began at 0000 UTC October 27, 1997 and ended at 1900 
UTC October 29. The simulation for the 30 km domain was started at 0000 UTC October 28 
and also ended at 1900 UTC October 29. Model data for the first 12 hours of each 
simulation were not utilized to allow for a sufficient “spin up ’ period. 
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Trajectories were calculated from the MM5 data using the kinematic method. That 
is, air parcels were advected using only the three dimensional wind field, without employing 
the isentropic assumption. The wind data at 26 sigma levels were interpolated to 24 1 equally 
spaced isobaric levels using a cubic spline procedure. Linear interpolation provided values 
between levels as well as at 5-minute time steps. Fuelberg et al. [1996] give further details 
about the trajectory model as well as a comparison between the kinematic and isentropic 
methods. 

Limitations of trajectories include incorrect placement of meteorological features by 
the input analyses. This is of particular concern over oceanic regions, such as the SONEX 
domain, where meteorological data are relatively sparse. However, numerical models such 
as the MM5 can generate the features that may be lacking in the coarser input data. 
Trajectories also are subject to the resolution of the wind data, since features of sufficiently 
small time or space scales will not be resolved [Doty and Perkey, 1992]. Resolution in our 
two domains was 1 hour and either 30 or 90 km. Thus, individual convective elements (e.g., 
their updrafts) still are not resolved. However, the effects of the convection are depicted at 
the scale of the grid. Finally, numerical limitations of the kinematic model itself affect 
overall trajectory accuracy [Stohl et al., 1995]. The kinematic method has been widely used 
in many meteorological/chemical applications [Draxler, 1991, Doty and Perkey , 1993, 
Krishnamurti et al., 1993, Chatfield et al., 1996, Garstang et al., 1996, and Bieberbach et al., 
this issue]. 

Arrival locations for the backward trajectories were at 1 -minute intervals along the 
flight path. Trajectory arrival times were the hour closest to the actual times of the in situ 
chemical observations. The trajectories were derived from a combination of the 30 and 90 
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km MM5 wind data. Specifically, the first 19-29 hours from the flight track were calculated 
from the 30 km (innermost) data grid, depending on time of arrival at the boundary of the 30 
km domain. The earlier 24-34 hours were calculated using the 90 km data. 

Locations of convection were determined using cloud-to-ground lightning data from 
the National Lightning Detection Network (NLDN) and Long Range Field (LRF) [ Cummins 
et al., 1998]. GOES-8 infrared imagery also was employed. 

2.2 Chemical Data 

The chemical data used in this study were taken from a 1 minute merged data set 
prepared at Harvard University. This data set contains 1-minute averages of all chemical, 
meteorological, and position data recorded during the flight, as well as model-derived values 
of certain chemical species. Singh et al. [this issue] describe the chemical measurements in 
detail. 

3. Meteorological Conditions 

High amplitude flow patterns dominated the SONEX period, including the day of 
Flight 10 [Fuelberg et al ., this issue]. Figure 3 contains winds and potential vorticity at 250 
hPa (-10.5 km) on 1200 UTC October 29 derived from the 1° x 1° ECMWF data (panels a,b) 
and the 90 km and 30 km MM5 simulations (panels c,d). 

Circulation patterns from the MM5 at 90 km resolution (Figure 3c) agree closely 
with those from ECMWF (Figure 3a), particularly the trough/ridge system near the NAFC. 
The cut-off low near the Azores also is similarly placed, although wind speeds from the 
MM5 are slightly greater (-10 m s' 1 ) than those in the ECMWF analyses. The departure site 
of Flight 10 (the Azores) was nearly coincident with this cyclone that had been drifting 
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slowly eastward after becoming detached from the main flow. The long wave trough and 
ridge axes (dashed lines) located south and west of the low, respectively, also are depicted 
similarly. The relatively deep trough south of Newfoundland was oriented over the 
southeastern United States two days prior to the flight. This system and its associated frontal 
boundaries produced extensive lightning from the central Gulf of Mexico to the coastal 
waters off the eastern United States. The two versions of the jet stream exhibit only slight 
differences in placement and wind speed. 

Potential vorticity (PV) (Figure 3b,d) provides a rigorous comparison between the 
ECMWF and MM5 data since PV requires derivatives of temperature and the horizontal 
wind components. The two versions of PV at 250 hPa exhibit close agreement, i.e., major 
features are similarly placed with only minor differences in magnitude. The MM5 generates 
an 1 1 PVU (1 PVU = lxlO' 5 K mb' 1 s' 1 ) maximum near the Azores at 30 km resolution 
compared to the 12 PVU maximum indicated by the ECMWF analysis. The pattern of 
isopleths indicates that the MM5 better resolves this tropopause undulation. The two PV 
maxima to the northeast and northwest of the main feature also are comparably represented 
in terms of magnitude, with maximum variations of only i PVU. 

Although only two parameters at 250 hPa are described here, we have compared the 
MM5 and ECMWF versions of many meteorological variables at various altitudes. Close 
agreement is found in all cases. Therefore, we believe that the MM5 simulation can be used 
confidently to understand the chemical signatures observed during Flight 10. Unless 
otherwise noted, future meteorological analyses will be derived from the MM5. 
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4. Results 

Results show that the complex chemical signatures observed during Flight 10 can be 
directly linked to the meteorological features described in Section 3. The stratospheric 
influence during the flight was among the greatest of the entire SONEX mission [ Thompson 
et al, this issue a]. The convective influence also was pronounced [Chen et al, this issue, 
Davis et al., this issue, Snow et al., this issue, and Thompson et al., this issue a], and a 
moderate amount of continental pollution was encountered. 

Three distinct chemical signatures were observed during Flight 10, and correlations 
between NO y and CO and between NO y and CN are useful in differentiating between them 
e.g., Thompson et al. [this issue, a], Jaegle et al [1996]. Figure 4 shows scatter diagrams of 
these three species during the flight. Stratospheric air is easily distinguished from lightning 
influenced or continental air since NO y is anti-correlated with both CO and CN in the 
stratosphere. Large NO y associated with large CO and large CN indicates lightning. 
Conversely, moderate NO y with moderate CO and CN denotes continental influence. 
Stratospheric and tropospheric regimes can be distinguished by the probability distribution of 
ozone mixing ratio (Figure 5). Although most ozone measurements during Flight 10 are less 
than 100 ppbv, numerous samples exceed the customary 100 ppbv stratospheric threshold, 
with some approaching 450 ppbv. 

Figures 6 and 7 are time series of selected chemical species measured along the 
path of Flight 10. Chemical signatures are denoted along the tops of the figures — two 
encounters of stratospheric air, three lightning signatures, one signature representing 
continental pollution, and a broad transition region along two flight legs. The bottom plot of 
each figure depicts flight altitude. Plate 1 shows the altitude and latitude of the flight and the 
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chemical signatures. The following sections describe chemical characteristics and transport 
paths associated with each signature. The two stratospheric portions of the flight (STR) are 
examined first. Next, we explore the three signatures associated with lightning (LT1, LT2, 
and LT3) and the pollution signal (POL). Finally, the transition region (TRS) is examined. 

4.1 Stratospheric Signature 

Well defined stratospheric air is encountered during Flight 10. Stratospheric 
penetration occurs on both the outbound and inbound legs when the DC-8 passes through the 
cutoff low near the Azores (Figure 3). During the outbound penetration, the stratospheric 
signal occurs between 1 128-1214 UTC (Plate 1, Figures 6,7), i.e., after the aircraft ascends to 

10. 1 km (267 hPa). Values of 0 3 and NO y increase sharply to -250 ppbv and 1000 pptv, 
respectively. Conversely, values of CO (< 40 ppbv), ultrafine CN (< 1000 cm' 3 ), and water 
vapor (< 100 ppmv) show dramatic decreases, while NO (< 100 pptv) and PAN (-40 pptv) 
decrease to a small fraction of NO y . Penetration during the return leg occurs at a higher 
altitude than before (11.9 km (202 hPa) between 1 725- 1 830 UTC), producing an even greater 
value of ozone (> 400 ppbv). 

The tropopause is relatively low near the cutoff low. This depression is depicted in 
the cross section of PV (Figure 8). If we assume that the stratospheric threshold is 3.0 PVU 
[Fuelberg et a/., this issue], the lowest altitude of the stratosphere is -7 km. The Dial 
(Differential absorption lidar) [Browell, 1989, 1991] ozone image (Plate 2) shows the 
downward protrusion of ozone rich stratospheric air. Specifically, 0 3 > 100 ppbv extends to 
-6 km during the return portion of the flight (right side of figure), with enhanced ozone (> 80 
ppbv) reaching nearly 3 km. Backward trajectories arriving at the stratospheric portion of the 
flight (not shown) remain at stratospheric altitudes throughout their 2-day histories. 
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4.2 Lightning Signatures 

Three regions of enhanced NO of varying intensity and duration are attributed to 
li ghtnin g. Trajectory analyses will show that each of these regions is associated with an 
extensive area of lightning that occurred along the East Coast of the United States 1 to 2 days 
prior to the flight. 

The longest chemical signature attributed to lightning (LT1) occurs from 1404 to 
1452 UTC (Figures 6,7). One should note that the bounds of LT1 encompass the entire flight 
leg at 9.5 km (292 hPa, Plate 1). The chemical signature of LT1 includes enhanced mixing 
ratios of NO (300-500 pptv), NO y (500-1000 pptv), and water vapor (200-300 ppmv), as well 
as enhanced ultrafine CN (10,000-17,000 cm' 3 ). However, CO is slightly reduced (~80 
ppbv). 

Figure 9 shows a horizontal plot of backward trajectories arriving at locations 
comprising LT1 . The trajectories pass over the Florida peninsula approximately two days 
prior to arrival, travel offshore of the Carolinas as they head northeast toward Newfoundland, 
and then turn southeast toward the flight track. To determine the extent to which parcels 
comprising LT1 were exposed to lightning, we superimpose trajectory and NLDN cloud-to- 
ground lightning data over 3 -hour intervals. Although intracloud lightning also produces 
nitrogen oxides, no data about its occurrence was available at the needed resolution. The 
trajectories experience two major encounters with cloud-to-ground lightning. The first 
occurs between 1200-1800 UTC October 27 (Plate 3) as they pass through the eastern Gulf 
of Mexico into Florida. The greatest exposure occurs between 1200-1500 UTC when the 
trajectories are nearly coincident with a concentrated area of cloud-to-ground flashes. 

Heights of individual trajectories during this encounter (47-51 hours back in Figure 10) vary 
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considerably, with the highest trajectories near 9.8 km and the lowest near 6.0 km. Satellite 
imagery (not shown) indicates cloud top temperatures in this area that correspond to heights 
of — 9- 1 2 km (~300 hPa to 200 hPa). Pickering et al. [1996] found that maximum convective 
outflow occurred near the tops of clouds that were investigated over Brazil. 

The trajectories comprising LT1 experience a second major encounter with 
lightning approximately 18 hours after the first (Plate 3). This occurs between 0600-0900 
UTC October 28 as the trajectories exit the northern edge of the flash region. Maximum 
cloud heights in this area (29-33 hours back in Figure 10) are ~12 km (200 hPa), while 
trajectory heights have risen to —7.0- 1 1.5 km (213-473 hPa). 

It is informative to examine lightning signature LT2 in the context of the LT1 
signature just described. LT2 is brief, extending only from 1322-1331 UTC (Figures 6,7). It 
occurs while the DC-8 ascends (Plate 1) from a leg at 8.8 km (319 hPa) to a leg at 10.7 km 
(243 hPa). Specifically, LT2 begins when the DC-8 ascends through the altitude of LT1 (9.5 
km, 292 hPa) in a location directly above LT1 . Values of NO, NO y , CN, CO, and water 
vapor for LT2 are very similar to those observed during LT1. In addition, the trajectories 
arriving at LT2 (Plate 4) have lightning encounters that are very similar in both time and 
space to those of the LT1 trajectories (Plate 3). Thus, LT2 appears to be a quick sampling of 
the same lightning event that produces LT1 . 

The LT2 signature ends part way through the 10.7 km leg at -1330 UTC (Figures 
6,7, Plate 1). This corresponds to an increase in 03 to ~80 ppbv, and a decrease in H2O to 
values similar to those encountered in the stratosphere. The DIAL image (Plate 2) indicates a 
downward protrusion of relatively large O3 in the region, and the horizontal analysis of PV 
(Figure 3b, d) shows that this location (~23 N, 35 W) lies along a narrow band of enhanced 



PV. A time series of flight level PV (Figure 9 in Fuelberg et al. [this issue]) indicates that 
values in this area increase to ~2.0 PVU. These factors suggest that the region exhibits a 
weak stratospheric influence. The only counter indication is that CO does not show the 
decrease that is expected of stratospheric air. 

The third lightning episode (LT3) occurs between 1646-1716 UTC as the DC-8 
heads back to the Azores after the stack maneuver. The signature begins at an altitude of 
10.3 km (245 hPa) while the aircraft is ascending from the 5.8 km (486 hPa) leg to the 1 1.9 
km (202 hPa) leg (Plate 1). The signal continues along the 1 1.9 leg until the second 
stratospheric penetration occurs (Section 4.1). The chemical signature of LT3 is the most 
intense of the three li ghtnin g events (Figures 6,7). LT3 exhibits the greatest NO (1500 pptv) 
and NO y (2000 pptv) of the entire flight. Values of ultrafine CN (30,000 cm' 3 ) also are 
greatest in this segment, while CO is relatively small ( — 80 ppbv). 

Trajectories arriving along the level leg of LT3 at 1 1.9 km are shown in Figure 1 1 . 
Their locations are not greatly different from those comprising LT1 (Figure 9). However, 
even these slight differences cause the LT3 trajectories to coincide more closely with the 
cloud-to-ground lightning locations, and for a longer period of time, than do those of the 
other lightning signatures. These lightning encounters are shown in Plate 5. The event 
begins near 1500 UTC October 27 as the first trajectories arrive at the southwestern edge of 
an extensive line of intense convection over the Gulf of Mexico. This encounter region, 
immediately north of the Yucatan Peninsula, contains the highest cloud tops of the entire line 
(> 12 km, < 200 hPa). At this initial time, only the southernmost trajectories coincide with 
the region of cloud-to-ground lightning. However, as the trajectories move eastward, their 
horizontal separation decreases, leading to a greater number of trajectory/lightning 
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encounters. Thus, between 1800-2100 UTC (Plate 5), a significantly greater number of 
trajectories are co-located with lightning. This first encounter ends at -2100 UTC October 
27, and over the next 12 hours (not shown), the number of lightning flashes decreases 
greatly. The trajectories pass over South Florida and emerge over the coastal waters of the 
Atlantic during this period. 

Cloud-to-ground lightning flashes increase dramatically over the western Atlantic 
Ocean after 0600 UTC October 28 (not shown), and the eastward moving trajectories pass 
through the southern part of this region between 0900-1500 UTC October 28 (Plate 5). One 
should note that almost all of the trajectories pass through the area of intense lightning at 
some time during this period. Over this six hour span, trajectory heights remain virt ual ly 
constant at near 1 2 km (200 hPa, not shown), corresponding approximately to cloud top 
level. 

To summarize, the strength of the LT3 chemical signature (Figures 6,7) apparently 
is due to the long period over which its trajectories are exposed to extensive lightning. This 
exposure for LT3 is greater than observed with the two other weaker li gh t nin g signatures 
(LT1-2). 

It should be noted that the high-resolution MM5 wind data were valuable for 
computing the trajectories arriving at LT3. For example, trajectories calculated from 
ECMWF data, where resolution (6-hourly, 1 degree) is coarser than the MM5, were not 
coincident with the observed lightning flashes at any time during their 2-day histories (not 
shown). 

Finally, a brief chemical "spike" occurs at -1220 UTC at the southern end of the 
10.1 km flight leg (Figures 6,7, Plate 1). It is located just south of the stratospheric region, 



near 32°N. The spike is characterized by slightly enhanced values of NO (-250 pptv), NO y 
(-600 pptv), and ultrafme CN (-8000 cm' 3 ). This signature also is consistent with a brief 
lightning encounter. Its altitude is between those of LT1-3. Thus, results indicated that the 
layer between -9-12 km located south of the stratospheric region is characterized by a 
lightning signature. 

43 Pollution Signature 

The pollution signature (POL) begins when the DC-8 descends to its lowest altitude 
of the flight (5.8 km, 486 hPa) (Figures 6,7, Plate 1). It occurs during the entire length of this 
leg, from -1528-1646 UTC (20°-28°N). Several species exhibit moderately enhanced values, 
including PAN (70-120 pptv), H 2 0 (500-1500 ppmv), and 0 3 (80 ppbv) in Figures 6,7 and 
C 2 H 2 (100-140 pptv), H 2 0 2 (250-450 pptv), and CH 3 OOH (75-150 pptv) (not shown). Values 
of NO and CN are not enhanced. 

Horizontal plots of the trajectories comprising POL are shown in Figure 12. 
Although their basic orientation is similar to those seen for previous signatures, the important 
difference is that these trajectories pass very near the heavily populated Northeast Corridor of 
the United States. Plate 6 focuses on the trajectories between 1500-2100 UTC October 27 
when they pass along the East Coast. One should note that there are virtually no encounters 
between the trajectories and cloud-to-ground lightning, i.e., the trajectories pass west of the 
lightning. This is consistent with the chemical signature. On the other hand, satellite 
imagery (Plate 6) indicates that the trajectories do pass near relatively deep clouds that do not 
contain cloud-to-ground lightning. Cloud locations at 2045 UTC are especially interesting 
(Plate 6c), i.e., a band of rain showers parallels the coast from North Carolina to Cape Cod. 
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Coldest cloud top temperatures range from -35° to -45°C, corresponding to heights of -7.5- 

9.5 km. The median trajectory altitude in this region is -8.6 km (between 36-42 hours back 
in Figure 13). Since this area is west of a cold front that is located offshore over the Atlantic, 
low level winds are from the northwest (not shown). Thus, it appears that surface based 
pollution from the Northeast Corridor is transported at low levels toward the line of 
convection, transported aloft by the convective updrafts, and then carried horizontally to the 
location of the POL signature. 

4.4 Transition Layer 

The layer between -7-9 km appears to be a transition region (TRS). It is sampled 
during two flight legs, i.e., first at 8.8 km (320 hPa) between 1245-1319 UTC, and again at 

7.6 km (376 hPa) between 1453-1526 UTC (Plate 1). A chemical “spike” is evident along 
the 8.8 km leg (indicated by an ellipse in Plate 1 and an arrow in Figures 6,7). The signature 
occurring at ~1310 UTC is characterized by the following values: NO (90 pptv), NO y (300 
pptv), and ultrafine CN (5000 cm' 3 ). Although these values are only slightly enhanced 
compared to the background, they are consistent with a weak lightning signature. 

Plate 7 shows that the trajectories comprising this signature experience only a slight 
encounter with cloud-to-ground lightning between 2100 UTC October 27 to 0300 UTC 
October 28 as they exit the coastal waters of the United States to the northeast. There are no 
close encounters during the periods that are not shown. These trajectories originate over 
Florida 49 hours earlier and then pass just south of Newfoundland before arriving at the 
flight track from the northwest (not shown). At the times when trajectory paths coincide with 



ligh tnin g locations, trajectory altitudes (not shown) are ~5. 7-6.2 km (-450-490 hPa), while 
satellite derived cloud top temperatures indicate heights near -7-9 km (-300-400 hPa). 

The 8.8 km flight leg containing the lightning “spike” is below the legs containing 
the major lightning signatures (LT1-3, Plate 1). And, trajectories arriving at the level of the 
spike pass over the coastal waters of the East Coast at lower altitudes than those arriving near 
LT1-3 (not shown). These lower levels appear to be below the main outflow of the extensive 
storms producing the intense cloud-to-ground lightning. Thus, the DC-8 may be sampling 
outflow from intermittent lower level lightning producing storms. This hypothesis is difficult 
to verify with the hourly satellite-derived cloud top data that are available. 

The lower altitude transition signature (TRS) begins at —1453 UTC when the DC-8 
descends to 7.6 km (376 hPa) (see Figures 6,7, Plate 1). The signature lasts -30 minutes 
(until 1526 UTC) when the DC-8 descends to its lowest leg containing POL. Plate 1 shows 
that this 7.6 km leg lies between the altitude legs containing LT1 (9.5 km) and POL (5.8 km). 
The chemical signature at 7.6 km is intermediate to those above and below. For example, 
note the values of NO (—100 pptv), NO y (-400 pptv), ultraflne CN (—8000 cm 3 ), PAN (—60 
pptv), water vapor (200-300 ppmv), and hydrogen peroxide (H 2 O 2 , -100 pptv, not shown). 

Horizontal locations of trajectories arriving along the 7.6 km TRS leg are shown in 
Figure 14. One should note that paths along the East Coast are intermediate to those of LT1 
(Figure 9), located directly above (at 9.5 km), and POL (Figure 12), located directly below 
(at 5.8 km). The LT1 trajectories are farthest offshore, encountering extensive cloud-to- 
ground lightning (Plate 3). Conversely, the POL trajectories are farthest west, avoiding 
virtually all of the lightning (Plate 6), but encountering continental pollution that was carried 
offshore by northwesterly winds behind a cold front. The intermediate TRS trajectories 
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experience some lightning encounters (not shown), but also appear to carry a pollution 
signature. The altitudes of the trajectories, when passing along the East Coast (Figure 15), 
support this hypothesis. One should note the wide range of trajectory altitudes (5.5-10 km, 

500-275 hPa) during passage from Cape Hatteras to east of Cape Cod between 36-42 hours 
prior to arriving at the DC-8. This contrasts with the more narrow range of altitudes 
exhibited by the LT1 and POL trajectories in this region (Figures 10, 13). This range of 
altitudes should encompass outflow from deep convection having extensive lightning as well 
as more shallow clouds with little or no lightning. 

The background chemical signature along the 8.8 km TRS leg near 1230 UTC (Figures 6,7, 
Plate 1) is not identical to that along the lower 7.6 km TRS leg. The upper TRS leg also contains the 
chemical spike that was attributed to a brief lightning encounter. These chemical differences are 
expected in a transition layer that is not fully mixed. 


5. Summary and Conclusions 

Flight 10 of NASA's Subsonic Assessment (SASS) Ozone and Nitrogen Oxide 
Experiment (SONEX) originated in Lajes, Azores (38.7°N, 27.1°W) and extended southwest 
to 19.9°N, 37. 1°W. A variety of chemical signatures were encountered during the flight. A 
stack maneuver along the southern half of the flight track provided samples at four different 
altitudes, i.e., 10.7, 9.5, 7.6, and 5.8 km. The horizontal length of these flight legs was 
relatively short, only -400-800 km. In this research we have examined the chemical 
signatures in detail, relating them to high resolution meteorological data from a mesoscale 


numerical model. 
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The National Center for Atmospheric Research (NCAR)/Penn State Fifth 
Generation Mesoscale Model (MM5) was used to create a high resolution three dimensional 
meteorological data set. The model was configured with two way interactive nesting, i.e., 
with a coarse (90 km) outer gnd and a finer (30 km) internal gnd. Each domain contained 26 
vertical levels. Analyses from the MM5 simulation agreed closely with global analyses. The 
MM5 data were available at hourly intervals and then used to create backward trajectories 
from the locations of the chemical signatures. We believe that the MM5-derived 
meteorological data provided a better representation of actual wind regimes and, 
consequently, more accurate trajectories than possible using global gndded analysis 
(typically 1 10 km, 6-hourly). This improved resolution is important since trajectories are 
very sensitive to slight errors in both wind direction and speed as well as the temporal 
frequency of the wind data. 

Four major categories of chemical signatures were discussed — stratospheric, 
lightning, continental pollution, and a transition layer. The strong stratospheric signal was 
encountered just south of the Azores. Values of Oj were greatly enhanced, while values of 
CO and water vapor were very small. The tropopause was relatively low in the region, 
extending down to ~6 km, due to a nearby closed low pressure circulation. 

Three chemical signatures were attributed to lightning. These signatures were 
located south of the lowered stratospheric region at relatively high altitudes (1 1.9, 10.7 and 
9.5 km). Two of the signatures occurred during the stack maneuver. Their chemical 
characteristics included enhanced values of NO, NO y , and water vapor, as well as enhanced 
ultrafine CN. Backward trajectories arriving at locations of the chemical signatures were 
related to locations of cloud-to-ground lightning. Results showed that the trajectories had 
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passed over regions of lightning located over the eastern Gulf of Mexico and off the 
southeast coast of the United States 1-2 days earlier. Two additional lightning signatures (the 
"spikes") were observed during legs at 8.8 and 10.1 km. 

The lowest leg of the stack maneuver (5.8 km) exhibited a chemical signature 
consistent with continental pollution. Species that were enhanced included PAN, water 
vapor, O3, C2H2, H2O2, and CH3OOH. Trajectories arriving at the pollution signature had 
passed over the highly populated Northeast Corridor of the United States. Since they passed 
to the west of the region of lightning that was farther offshore, no lightning signature was 
indicated. Surface based pollution apparently was lofted to the altitudes of the trajectories by 
convective clouds along the East Coast that did not contain lightning. 

Finally, a transition layer was detected along two flight legs between ~7-9 km. The 
chemical signature at 7.6 km was intermediate to signatures located at higher and lower flight 
legs. Species having intermediate values included NO, NO y , CN, PAN, water vapor, and 
hydrogen peroxide. The lightning signature was located directly over the 7.6 km leg at an 
altitude of 9.5 km, while the pollution signature was located directly below at 5.8 km. 
Trajectories arriving in this region had passed over locations off the East Coast that were 
between those of the lightning and pollution signatures. Thus, they probably were impacted 
by both lightning and pollution. These trajectories exhibited a wide range of altitudes when 
passing over the region of influence. A second leg of transition was located at 8.8 km. The 
chemical signature along this leg was not as pronounced as that at 7.6 km. 

In summary, the results show that high resolution numerical modeling can provide 
meteorological data that are useful for understanding complex chemical signatures. This 
capability was especially evident during the aircraft’s stack maneuver, when four different 
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chemical signatures were co-located at different altitudes, i.e., pollution, transition, lightning, 
and stratospherically influenced. 
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Plate Captions 


Plate 1. Altitude profile of Flight 10 as a function of latitude. Chemical signatures are 
indicated by colored segments as follows: Red-stratospheric, blue-lightning, green- transition, 
and orange-pollution. Lightning signatures within the transition leg and at the end of the 
stratospheric leg are circled. 

Plate 2. Cross section of DLAL-derived O 3 mixing ratio. The O 3 scale (ppbv) and time 
(UTC) are give at the top while latitude and longitude are given at the bottom. Altitude in 
kilometers is given along the sides. 

Plate 3. Lightning event 1 . Three hour trajectory paths are indicated by blue circles. 
NLDN cloud-to-ground lightning flashes are indicated in red. The time segment is given at 
the top of each panel. 

Plate 4. Lightning event 2. Three hour trajectory paths are indicated by blue circles. 
NLDN cloud-to-ground lightning flashes are indicated in red. The time segment is given at 
the top of each panel. 

Plate 5. Lightning event 3. Three hour trajectory paths are indicated by blue circles. 
NLDN cloud-to-ground lightning flashes are indicated in red. The time segment is given at 
the top of each panel. 
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Plate 6. a) GOES-8 infrared satellite image at 1445 UTC October 27, 1997. 

b) Three hour trajectory paths (blue) and NLDN cloud-to-ground lightning 
strikes (red) for the period of 1500-1800 UTC October 27, 1997. 

c) GOES-8 infrared satellite image at 2045 October 27, 1997. 

d) Three hour trajectory paths (blue) and NLDN cloud-to-ground lightning 
strikes (red) for the period of 1800-2100 UTC October 27, 1997. 

Plate 7. Transitional segment. Three hour trajectory paths are indicated by blue circles. 
NLDN cloud-to-ground lightning flashes are indicated in red. The time segment is given at 
the top of each panel. 


Figure Captions 

Figure 1. DC-8 flight track on October 29, 1997. 

Figure 2. Domains for the MM5 simulation. The outer perimeter indicates the boundary 
of the 90 km domain. The inner box represents the area of the 30 km domain. The track of 
Flight 10 is indicated within the 30 km domain. 

Figure 3. a) ECMWF 250 hPa streamlines and isotachs at 1200 UTC October 29, 1997. 
Wind speeds are contoured at intervals of 10 m s' 1 . 
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b) ECMWF 250 hPa potential vorticity at 1200 UTC October 29, 1997. Contours 
are 1 PVU. 

c) MM5 250 hPa streamlines and isotachs at 1200 UTC October 29, 1997. Wind 
speeds are contoured at intervals of 10 m s' 1 . 

d) MM5 250 hPa potential vorticity at 1200 UTC October 29, 1997. Contours are 1 
PVU. The axis of a cross section in Figure 8 is shown. 

Figure 4. a) Correlation plot of NO y (pptv) versus CO (ppbv). The three main signatures 
of Flight 10 are labeled. 

b) Correlation plot of NO y (pptv) versus the logarithm of ultrafine CN 
(cm' 3 ). The three main signatures of Flight 10 are labeled. 

Figure 5. Probability distribution function of O 3 for Flight 10. Ozone values of 100 ppbv 
or greater indicate stratospheric air. 

Figure 6 . Times series of chemical species and DC -8 altitude along Flight 10. Mixing 
ratios of NO (pptv), NO y (pptv), CO (ppbv), and O3 (ppbv) along the flight track are given in 
the first four panels. The bottom panel indicates the altitude of the DC- 8 . Chemical 
signatures are indicated on the top panel (STR-stratospheric, TRS-transition, LT-lightning, 
and POL-pollution). Lightning events within the transition zone and at the end of the 
stratospheric leg are indicated with arrows. 
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Figure 7. Times series of chemical species and DC-8 altitude along Flight 10. 
Concentrations of ultrafine CN (cm 3 ) and mixing ratios of PAN (pptv) and H 2 O (ppmv) 
along the flight track are given in the first three panels. The bottom panel indicates the 
altitude of the DC-8. Chemical signatures are indicated on the top panel (STR-stratospheric, 
TRS-transition, LT-lightning, and POL-pollution). Lightning events within the transition 
zone are indicated with arrows. 

Figure 8 . Cross section of MM5 derived potential vorticity (PVU) through the cut-off low 

pressure system at 1200 UTC October 29, 1997. Altitude (km) and horizontal distance are 
indicated. The axis is shown in Figure 3d. 

Figure 9. Backward trajectories arriving at LT1 (50-51 hours back). Arrival time at flight 
level is 1400-1500 UTC October 29, 1997. Arrows along trajectory paths indicate locations 
at six hour intervals. 

Figure 10. Trajectory altitudes (km) as a function of time (hours back) for arrivals at LT1 
(1404-1452 UTC). 

Figure 11. Backward trajectories arriving at LT3 (53 hours back). Arrival time at flight 
level is 1700 UTC October 29, 1997. Arrows along trajectory paths indicate locations at six 


hour intervals. 
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Figure 12. Backward trajectories arriving at POL (51-53 hours back). Arrival time at 
flight level is 1500-1700 UTC October 29, 1997. Arrows along trajectory paths indicate 
trajectory locations at six hour intervals. 

Figure 13. Trajectory altitudes (km) as a function of time (hours back) for arrivals at POL 
(1527-1646 UTC). 

Figure 14. Backward trajectories arriving at portion of the transition (TRS) segment (5 1 
hours back). Arrival time at flight level is 1500-1700 UTC October 29,1997. Arrows along 
trajectory paths indicate locations at six hour intervals. 

Figure 15. Trajectoiy altitudes (km) as a function of time (hours back) for arrivals at the 
1453-1526 UTC portion of TRS. 
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Plate 1. Altitude profile of Flight 10 as a function of latitude. Chemical signatures 

are indicated by colored segments as follows: Red-stratospheric, blue-lightning, green- 
transition, and orange-pollution. 
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Plate 2. Cross section of DIAL-derived 0 3 mixing ratio. The 0 3 scale (ppbv) and 

time (UTC) are give at the top while latitude and longitude are given at the bottom. 
Altitude in kilometers is given along the sides. 
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Plate 4. Lightning event 2. Three hour trajectory paths are indicated by blue 

circles. NLDN cloud-to-ground lightning flashes are indicated in red. The time segment 
is given at the top of each panel. 
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Plate 5. Lightning event 3. Three hour trajectory paths are indicated by blue 

circles. NLDN cloud-to-ground lightning flashes are indicated in red. The time segment 
is given at the top of each panel. 










Plate 6. a) GOES-8 infrared satellite image at 1445 UTC October 27, 1997. 

b) Three hour trajectory paths (blue) and NLDN cloud-to-ground lightning 
strikes (red) for the period of 1500-1800 UTC October 27, 1997. 

c) GOES-8 infrared satellite image at 2045 October 27, 1997. 

d) Three hour trajectory paths (blue) and NLDN cloud-to-ground lightning 
strikes (red) for the period of 1 800-2 1 00 UTC October 27, 1997. 




971 027 21 00-971 028 0000 UTC 



Plate 7. Transitional segment. Three hour trajectory paths are indicated by blue 
circles. NLDN cloud-to-ground lightning flashes are indicated in red. The time segment 
is given at the top of each panel. 
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Figure 1. 


DC-8 flight track on October 29, 1997. 




Figure 2. Domains for the MM5 simulation. The outer perimeter indicates the 
boundary of the 90 km domain. The inner box represents the area of the 30 km domain. 
The track of Flight 10 is indicated within the 30 Ion domain. 





Figure 3. a) ECMWF 250 hPa streamlines and isotachs at 1200 UTC October 29, 
1997. Wind speeds are contoured at intervals of 10 m s' 1 . 

b) ECMWF 250 hPa potential vorticity at 1200 UTC October 29, 1997. 
Contours are 1 PVU. 

c) MM5 250 hPa streamlines and isotachs at 1200 UTC October 29, 1997. 
Wind speeds are contoured at intervals of 10 m s' 1 . 

d) MM5 250 hPa potential vorticity at 1200 UTC October 29, 1997. 
Contours are 1 PVU. The axis of a cross section in Figure 8 is shown. 
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Figure 4. a) Correlation plot of NO y (pptv) versus CO (ppbv). The three main 
signatures of Flight 10 are labeled. 

b) Correlation plot of NO y (pptv) versus the logarithm of ultrafme CN 
(cm' 3 ). The three main signatures of Flight 10 are labeled. 
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Figure 6. Time series of chemical species and DC-8 altitude along Flight 10. 
Mixing ratios of NO (pptv), NO y (pptv), CO (ppbv), and Oj (ppbv) along the flight track 
are given in the first four panels. The bottom panel indicates the altitude of the DC-8. 
Chemical signatures are indicated on the top panel (STR-stratospheric, TRS-transition, 
LT-lightning, and POL-pollution). Lightning events within the transition zone and at the 
end of the stratospheric leg are indicated with arrows. 
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Figure 8. Cross section of MM5 derived potential vorticity (PVU) through the cut- 
off low pressure system at 1200 UTC October 29, 1997. Altitude (km) and horizontal 
distance are indicated. The axis is shown in Figure 3d. 
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Figure 1 1 . Backward trajectories arriving at LT3 (53 hours back). Arrival time at 
flight level is 1700 UTC October 29, 1997. Arrows along trajectory paths indicate 
locations at six hour intervals. 






Figure 12. Backward trajectories arriving at POL (51-53 hours back). Arrival time at 
flight level is 1500-1700 UTC October 29, 1997. Arrows along trajectory paths indicate 
trajectory locations at six hour intervals. 
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Ficure 14 Backward trajectories arriving at portion of the transition (TRS) segment 
(Slhoursback). Arrival time at flight level is 1500-1700 UTC October 29,1997. Arrows 
along trajectory paths indicate locations at six hour intervals. 










